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ABSTRACT 

IRAS 20050+2720 is young star forming region at a distance of 700 pc without apparent high 
mass stars. We present results of our multiwavelength study of IRAS 20050+2720 which includes 
observations by Chandra and Spitzer, and 2MASS and UBVRI photometry. In total, about 300 YSOs 
in different evolutionary stages are found. We characterize the distribution of young stellar objects 
(YSOs) in this region using a minimum spanning tree (MST) analysis. We newly identify a second 
cluster core, which consists mostly of class II objects, about 10' from the center of the cloud. YSOs 
of earlier evolutionary stages are more clustered than more evolved objects. The X-ray luminosity 
function (XLF) of IRAS 20050+2720 is roughly lognormal, but steeper than the XLF of the more 
massive Orion nebula complex. IRAS 20050+2720 shows a lower Nh/Ak ratio compared with the 
diffuse ISM. 

Subject headings: circumstellar matter - infrared: stars - Stars: formation - Stars: pre-main sequence 
- X-rays: stars 

1. INTRODUCTION 

Star formation occurs in dense and cool molecular 
clouds, which collapse under their own gravity. We 
observe those clouds over a range of masses, but the 
majority of stars forms in clusters with more than a 
hundred members (|Carpenterl 120001 : lLada fc Ladall2003t 
iPorras et~a l. 2003). Thus, studying star formation in 
those clusters is an important step to understand the 
history of the stars we currently observe, and the stellar 
population of the outer gala xy as a whole (for g eneral 
reviews of star formation see lReipurth et aT1l2007t ). 

In the early stages the proto-star is hidden by its par- 
ent circumstellar material and thus can only be observed 
as far-IR emission from the warm, dusty envelope (class 
0). This envelope becomes less dense in class I objects. 
They often drive powerful outflows and carve holes in 
their envelopes so radiation can escape. Eventually, the 
envelope disperses and the star becomes visible as a clas- 
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sical T Ta uri star (CT TS) or "class II" source in IR clas- 
sification (|Ladalll987l ). The disk still causes an IR excess 
over photospheric stellar emission. On the star itself ac- 
cretio n shocks and coronal activity lead to X-ray emis- 
sion tiLamzinl 119981 iKastner et al.l 120021 iGiinther et~all 
[20071 IGiintherl 1201 II) . Lat er, the IR excess va nishes as 
the gas disk is dissipated (jlnglebv et al.l [20121 ) . In this 
stage (weak-line T Tauri stars = WTTS) cluster mem- 
bers cannot be distinguished from main-sequence (MS) 
stars by IR observations. One method to identify them is 
through their X-ray luminosity, which is far higher than 
for most older stars (for a review of X-ray properties 
of young stellar objects (YSOs) before high-resolution 
X-ray spectroscopy with Cha ndra or XMM-Newton see 
iFeigelson fc Montmerrelll999D . 

Many details of the star formation process in clusters 
are still under debate despite significant observational 
progress. For example, it is unclear if in an undisturbed 
cluster high-mass or low-mass stars form first. However, 
the theoretical assumption of a cluster evolving with min- 
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imal influence from the rest of the gal axy may be far 
from realistic. For example, in p Oph Ide Geusl (|1992ft 
suggested that the expanding shell of the upper Sco star 
forming region travelled through the cloud and triggered 
the star formation. Events like that would cause a non- 
uniform age distribution within a cluster. 

In this paper we study the young cluster 
IRAS 20050+2720, which is lo cated at a distanc e 
of 700 pc based on CO velocities (jWilking et all 11989ft . 
It is found at a galactic longitude where the Cygnus Rift 
and the Cygnus X region are distinguishable, therefore, 
there is no distance ambiguity. We use this distance 
throughout the article. There appear to be no massive 
stars in IRAS 20050+2720, thus the intensity of the 
ambient radiation should be small and we can study the 
evolution of late-type stars in the absence of external 
irradiation on their disks. 

IRAS 20050+2 720 is associated with a water maser 
(iPalla et al.lll99il ) and a methanol maser (jFontani et alJ 
12010ft. Molecular lin e emission indicates mass infall 
(|Gregersen et al.lll997tlChoi et al.lll999f ). Radio observa- 
tions reveal a complex structure of several lobes close to 
the nominal source position of IRAS 20050+2720, which 
has been identified as two or three jets from different 
proto-stellar sources (iBachiller et al.lll995l ; ICodella et al.l 
Il999t iBeltran et all 12008ft . The luminosity of this re- 
gion is estimated a s 388 L m from th e IRAS fluxes 
(|Molinari efalll 1996ft . IChini et al.l (|200lD obtained mm- 
maps of the regions, estimatinge a total gas mass of 
200 M Q within 65" (0.2 pc) radius. 

The first attempt to establish an IR classification of in- 
dividual stars in IRAS 20050+2720 identified about 100 
clust er members, abo ut half of which show an IR ex- 
cess (jChen et al.lll997l ). While they estimate an average 
cluster age of 1 Myr, they interpret the radio lobes and 
the strong reddening as signatures for a more recent star 
forma tion event in the last 0.1 Myr. iGutermuth et all 
(2005) confirmed these ideas. 

With Spitzer it is possible to get very accurate pho- 
tometry of star forming regions and classify individual 
objects according to their IR spe ctral energy distribution 
(SED). IGutermuth et"aH (12009ft identified 177 YSOs in 
IRAS 20050+2720, wh ich they grouped i n two distinct 
cores. Compared with IChen et all (|1997l ) they identify 
more absorbed cluster members and show that the earlier 
definition of subclusters is partially due to an observa- 
tional bias. In this paper, we extend the spatial coverage 
with an additional Spitzer observation and add informa- 
tion from a Chandra observation and, surprisingly, the 
first dedicated optical photometry. 

Our goal is to compare the clustering and the X-ray 
properties of IRAS 20050+2720 with other star forming 
regions and to search for internal gradients of density and 
age that could indicate triggered star formation. 

In Sect. [2] we describe the observations of 
IRAS 20050+2720. In Sect. [3] we identify disk-bearing 
YSOs by the IR colors and disk-less YSOs based on 
their X-ray emission. Our results in Sect. |4] describe the 
spatial distribution of YSOs and the optical photometry 
allows us to place individual members on evolutionary 
tracks. We show the X-ray properties and discuss the 
Ak/Nh ratio. We end with a short summary in Sect. [5l 

2. OBSERVATIONS 



TABLE 1 
Log of observations 



Facility Id instrument exposure date 



In this section we describe the datasets we use for our 
analysis. They are summarized in table [TJ 

2.1. Spitzer 

Spitzer has observed this region several times. We use 
data from two initial observa tions performed in 2004 , 
which are described in detail in IGutermuth et all (2009). 
This dataset misses an extended population beyond the 
cluster core, so follow-up observations were performed 
in 2007 as part of the Spitzer "Clusters Near Clusters" 
GTO program (Program ID 40147) (jMasiunas et all 
12012ft . In 2004 t he cluster ce nter was imaged with IRAC 
(3.6 - 8.0 ^m, iFazio et al.l I2004T ) with a 3 x 4 raster 
map of four dithered images per position in the high 
dynamic range mode. In 2007 two additional 3x4 
raster maps were obtained to the north-east and south- 
west of the original pointing, which overlap partially 
with the initial observation. Table [T] gives the expo- 
sure time per pixel for each observati on. MIPS mo- 
saics (24 — 160 fim, iRieke et all 12004ft were taken at 
medium scan rate with full width scan stepping. The 
observation in 2004 contains three scanning stripes, one 
stripe towards the south-east was added in 2007. Be- 
cause of the low angular resolution the 70 and 160 /im 
data cannot be used for individual sources in a dense 
region like IRAS 20050+ 2720. All dataset s are r educed 
and combined following IGutermuth et all (|2009ft using 
the Clust erGrinder IDL packag e, which is described in 
detail bv lGutermuth et aTT (|2b08). In short, the standard 
basic calibrated data products (BCD) from the Spitzer 
Science center pipeline are treated for bright source ar- 
tifacts and combined into a mosaic image, eliminating 
outliers to remove cosmic rays. The mosaic of individual 
AORs are shifted to match the 2MASS astrometry. Au- 
tomated source detection and aperture photometry are 
carried out on the mosaics. 

Compared wi t h th e reduction described by 
IGutermuth et al.l (|2009l ) our new data reduction is 
improved in several ways to in crease the photometric 
accuracy (M asiunas et alJ 12012ft . Our data are now 
processed with BCD S18.5 which contains a muxbleed 
correction, resulting in a cleaner mosaic close to the 
bright sources in the center of the cluster. Also, we 
now use a pixel scale of 0.86", which is 1/V2 of the 
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native pixel size. This finer grid gives better centroid 
positions; it also allows to better separate close sources. 
Sect. 13.31 describes how the improved reductions effect 
the detection efficiency. 

2.2. Chandra 

IRAS 20050+2720 was observed with Chandra using 
the VFAINT mode of the ACIS-I detector. As indi- 
cated in table[T]the cluster was observed in three different 
epochs. The position angle varies between observations, 
thus the coverage of the field is non-uniform on the edges. 
The individual exposures are merged an d processed with 
the p ipeline of the ANCHORS project ([Spitzbart et al.l 
|2005|) . Source detecti on is performed with th e wavdetect 
algorithm in CI AO (jFruscione et al.1 120061 ). We keep 
sources with a significance above 2a. For each source 
with more than 20 counts, event lists and response files 
are generated. We fit a single Ray mond-Smith thermal 
emission model ([Foster et alJl2010f i. where the absorbing 
column density Nh, the temperature T and the volume 
emission meassure are free paramet ers. The metalicity 
is fix ed at 0.3 compared with solar ([Anders fc Grevessa 
1989). Also, for eac h source a lightcurve is gen erated and 
its Gregory-Loredo ([Gregory fc Loredolll992D statistic is 
calculated to characterize the source variability but this 
will be discussed in a later paper. The spectral proper- 
ties presented here are effectively time-averaged over all 
three observations. 

In total 368 sources are detected with Chandra. Ta- 
ble [2] lists the position of the X-ray sources and their 
properties (fx, Nh and kT and observed flux in the en- 
ergy range 0.3-8.0 keV) with statistical errors. Also, we 
calculated the unabsorbed luminosities assuming a dis- 
tance of 700 pc. In addition, this table contains the ID 
number of the source in the merged source table (table[3]) 
and the classification (see Sect. I3.2ll3.4[) . Table [5] is or- 
dered by the ID number from the merged source table, 
which in turn is ordered by the classification. 

2.3. Optical photometry 

Optical photometry of IRAS 20050+2720 was obtained 
at the Fred Lawrence Whipple Observatory (FLWO) 
with the 1.2 m telescope and the Keplercam instrument. 
Keplercam is a camera with a Fairchild "CCD 486" de- 
tector. Observations were taken with staggered exposure 
times of 6, 60 and 300 s, so that photometry for brighter 
stars can be performed on frames with shorter exposure 
times. Keplercam is equipped with filters for different 
photometric systems, here we use observations taken in 
u 1 B hVht' Y '. The Bjj and Vh filters belong to the Harris 
set, the u', r' and i' filters are part of a Sl oan Digital Sky 
Survey (SPSS) set dFukueita et al.lfl996h . 

We performed the data reduc tion in Pyraf, a Python 
interface to IRAF ([Todvl I1993 D Photometr y is done 
with PSF-fitting in DAOPHOT (|Stetsonlll987l) as imple- 
mented in IRAF. Data and flatficld are bias-corrected, 
then the data images are divided by the flatfield. 

The absolute pointing accuracy of the FLWO telescope 
is insufficient to match source positions with the other 
datasets directly. Therefore, we first correct the WCS of 
each individual exposure with the IRAF task msccmatch, 
which minimizes the residuals b etween sources in ou r 
images and the 2MASS catalog ([Skrutskie et al.ll2006D . 



Sources are extracted independently for different expo- 
sures and source lists are merged with a matching tol- 
erance of 3", which is of the order of the remaining as- 
trometric errors. Properties of bright sources are taken 
from shorter exposures to avoid the non-linear regime of 
the detector. 

To perform absolute photometry extinction and zero 
point of the instrumental magnitu de are calibra ted with 
observations of Landolt field 01 dStetsonl 120001 ). which 
was reduced analogues to the data of IRAS 20050+2720 
itself. We fit the instrumental magnitudes in v! BnVHr'i' 
to the catalog values for UBVRI assuming an offset, 
a linear color term and a linear extinction term. Ten 
sources in the standard field are extremely red (1.15 > 
R — I > 2.0), but we do no t find the break in t he trans- 
formation law indicated in iSmith et al.l (|2002f) for very 
red sources. Instead the full sample fits the relation 
R — I oc r' — i' within the observational uncertainties. 
We remove the faintest stars from our final sourcelist, so 
that the typical error on the flux is about 10%. 

The optical photometry is supplemented by data from 
IPHAS, the Isa ac Newton Telesc ope (INT) Photomet- 
ric Ha Survey ([Drew et all 120051) . IPHAS is a survey 
of the northern Galactic Plane in broad-band r', i', and 
narrow-band Ha filters to a depth of r 1 ~ 20 (10cr). Using 
the I PHAS Initial Data Release ([Gonzalez-Solares et al.l 
2008) we extracted all stellar sources (morphological 
classification = -1) within 15' of the center of IRAS 
20050+2720. Because IPHAS images the sky twice (with 
pointings offset to cover the Wide Field Camera chip 
gaps) most sources are observed twice and so these were 
cross-matched and the mean magnitudes calculated. We 
remove all sources with photometric uncertainties in any 
band above 0.1 mag. This results in a total of ~1200 
sources. 

3. DATA ANALYSIS 
3.1. Source matching 

The different infrared observations are matched as in 
iGutermuth et al.l (|2008|) . Uncertainties on the coordi- 
nates are 0'.'25 (lc). We then add the X-ray data, match- 
ing X-ray sources and IR sources, if they are closer than 
the large axis of the 3<r error ellipse from wavdetect. If 
more then one IR source fullfills this criterion, we match 
the closest X-ray and IR-source (less than 5 cases). If 
two X-ray sources are matched to the same IR source we 
keep the match, which is closer to the IR centroid and 
list the remaining X-ray source as unmatched (1 case). 
Optical sources are identified with IR or X-ray sources 
using a constant matching radius of 3". 

The final list of cross-matched sources is published on- 
line as an electronic table, table [3] describes its content. 
This table also contains a classification (see Sect 13. 2113. 4p 
and an estimate of the total optical extinction Ay of the 
cloud at the sour c e pos ition from the extinction map of 
IGutermuth eTail ((20091) . This map is based on H — K 
colors of 2MASS source s, assuming an average intrinsic 
color of (H - K) = 0.2 (|Gutermuth et al.ll2005h . 

3.2. Class I and II sources 

Class I and II sources are surrounded by circumstel- 
lar material, thus they have a distinct IR signature. We 
use the same classification method as in lGutermuth et all 
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TABLE 2 

X-RAY SOURCES AND SPECTRAL FITS 



ID 




a 
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Note. — This table shows only X-ray sources, that are identified 
detected X-ray sources can be found in the electronic version. 

(|2009(1 . which is based on a series of cuts in color-color di- 
agrams (CCD) and color-magnitude diagrams (CMD) to 
remove various contaminants (galaxies with PAH emis- 
sion, AGN, Herbig-Haro objects). The classification is 
primarily based on IRAC fluxes; however, if either the 
5.8 ^m or the 8.0 /im flux is not available, then the 
2MASS bands can be used to estimate Ak and dered- 
den the remaining IRAC data to perform the classifica- 
tion. The deredden i ng is based on the extinction law 
from iFlahertv et al.l (|2007T > and the reddening loci from 
iMeyer et al.l ( 1997h . In a final step, sources which are 
highly reddened and thus invisible at short wavelengths, 
but bright in the MIPS 24 pirn band are added. Likely, 
they are deeply embedded protostars, but their evolu- 
tionary status cannot be confirmed without submillime- 
ter data. In the tables they are denoted as class I* 
to distinguish them from less embedded class I sources. 
Sources with less surrounding material but an IR ex- 
cess in all bands due to a disk are classified as class II; 
those where the IRAC bands follow the energy distribu- 
tion of a stellar photosphere, but an excess is found in 
MIPS 24 fxm are marked as class II*. These sources are 
sometimes called transition objects because an excess is 
present only at longer wavelength and therefore larger 
disk radii. In an evolutionary sequence this can be inter- 
preted as a transition phase between disk-bearing class II 
objects and MS stars. Sources compatible with stellar 



as YSOs and have valid spectral fits. The full table including all 

photospheres are marked as class III candidates. These 
could either be young pre-main sequence WTTS or main- 
sequence (MS) stars and are further classified based on 
optical and X-ray properties. Table |4] gives the total 
number of detected IR sources. The dataset contains 57 
class I sources, about 20% of them (10) are detected in 
X-rays, and 183 class II sources, 30% (53) of them are 
detected in X-rays. 

3.3. Improved Spitzer data reduction 

In Sect. 12.11 (see also iMasiunas et al.1 12012| ) we de- 
scrib ed several changes in ou r data reduction compared 
with iGutermuth et all (|2009t ). The improved photom- 
etry allows us to detect additional YSOs and in some 
cases objects are reclassified. The number of detected 
class I sources rises from 41 to 51 in the central part of 
the 2004 observations. This field is not covered by the 
new observations, thus, the higher number is entirely due 
to an improved data treatment. The number of class II 
sources rises from 95 to 105 and the number of deeply 
embedded sources doubles from 6 to 11. In both cases 2 
sources are classified as transition disks (class II*). 

3.4. Class III sources 

The main goal of this study is to extend the treat- 
ment of YSOs to include young disk-less cluster mem- 
bers. While class I and class II objects can be identified 
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TABLE 4 
Number of IR sources 
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a : WTTS and MS stars both show stellar photospheric colors in 
the IR. b : Foreground stars. 

by their IR excess, class III sources have lost their enve- 
lope and their disk. In the IR they are indistinguishable 
from MS-stars. However, because of their low age they 
are still fast rotators and consequently luminous X-ray 
emitters. Our set of X-ray sources which do coincide 
with an IR source consistent with stellar photospheric 
IR colors, thus provides a list of strong candidates. All 
IR sources consistent with stellar photospheric IR col- 
ors but without an X-ray counterpart are reclassified as 
(main-sequence) "star" in table [3] 

3.4.1. Optical data 

Different methods can be applied to identify fore- 
ground sources in the sample. T he three brighte s t star s 
are spectroscopically classified bv lNesterov et al.l (|1995| ). 
thus their distance is easy to determine. Beyond that we 
rely on the IPHAS catalog and our own photom etry and 
use a procedure similar to Wright et alj (|2010t ). Fig. Q] 
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Fig. 1. — Color-color diagram of all IPHAS sources within 15' of 
the center of the cluster. X-ray sources are marked w ith squares. 
Overplotted is an unreddened and a reddened MS from D rew et all 
(2005). Sec the electronic edition of the Journal for a color version 
of this figure. 



shows the IPHAS CCD for all sources within 15' of the 
cluster center. Overplotted are all X-ray detections. 

Class II stars and young class III objects are typically 
Ha emitters, thus they are found above the MS. Giants 
and AGB stars are also found in this region, the latter ex - 
tending towards much redder colors (| Wright et al.ll2008h . 
There are eleven sources with r' — Ha significantly above 
the MS (more than three times the measurement uncer- 
tainty, but at least 0.2 mag) in figured] We exclude the 
object at r' — Ha = 1.5 and r' — i' = 4 which probably is 
an AGB star (jWright et al.l l2008l . Four of the remaining 
ten sources are outside the FOV of IRAC bands 1 and 3 
(ID 23216, 25132, 25423, 27352 in table O, one is unde- 
tected in IRAC band 3 and 4 (ID 12295) and thus those 
sources do not have an IR classification. Four sources 
coincide with class II sources (ID 92, 100, 105, 146). The 
remaining source is IR-classified as a star (ID 1645). 

Sources which are reddened by dust are located be- 
low and to the right of the MS, thus we expect cluster 
members and background sources in this region. Fore- 
ground objects should cluster around the MS; however, 
reddened young stars with Ha in emission could also oc- 
cupy this region, e.g. a young star with Ay- = 6 and Ha 
emission such that r' — Ha is about 0.3 mag higher than 
for a MS star would be very close to the unreddened 
MS in figure [TJ In this step we take all X-ray sources 
within 0.2 mag in either color of the MS as potential 
foreground sources. Next we plot those sources on a J 
vs. J — H color-magnitude-diagram in Fig. [2] and com- 
pare their position with a MS at 700 pc. In this diagram 
reddening moves sources of intermediate mass roughly in 
parallel to the MS, thus the moderate reddening experi- 
enced by foreground stars does not shift them into the 
region where the YSOs of IRAS 20050+2720 are located. 
Given the spread of all IR sources in the IPHAS CCD, 
it seems very unlikely that a significant fraction of YSOs 
with Ha emission coincides by chance with the unred- 
dened MS in figure Q] and also with the MS on the J vs. 
J — H diagram. Thus we remove nine sources which are 
compatible with unreddened MS stars at distances be- 
tween 100 and 700 pc on both diagrams from our list of 
young class III candidates. 
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Fig. 2. — J vs. J-H color-magnitude diagram. Filled colored 
symbols show X-ray sources with known parallax (green pentagon), 
those which are close to the MS in Fig. [I] and thus are potential 
foreground sources (red circl es) and oth e r X-r ay sources (yellow 
squares). The early MS from Siess ct al. (2000) is ovcrplotted for 
a distance of 700 pc. See the electronic edition of the Journal for 
a color version of this figure. 
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Fig. 3. — Distribution of X-ray sources. The contours show 
the extinction from a background extinction map. The countour 
levels are Ay = 2, 6, 10, 14 mag a fter subtracting a background 
of Ay = 2 (Gutcrmuth ct al. 2009). The gray scale image shows 
IRAC 8 fim emission. See the electronic edition of the Journal for 
a color version of this figure. 

3.4.2. Other X-ray sources 

MS stars in front of the cloud are detected in the opti- 
cal and IR down to spectral type M. Thus, the remain- 
ing unidentified sources are mostly background objects, 
which should be randomly distributed over the Chandra 
field-of-view. Fewer X-ray sources are expected towards 
the edges because the effective area of the ACIS detec- 
tor is lower here than in the center. Additionally, we 
expect fewer background X-ray sources where the clus- 
ter is densest, because the extinction hides weak back- 
ground sources. Figure [3] shows the distribution of all 
X-ray sources. 

Squares mark all sources which have been classified 
above, i.e. they belong to one of the following groups: I*, 
I, II, II*, young class III or MS. Small circles and triangles 
mark all as yet unidentified sources. They could be back- 
ground sources like AGN, high-mass X-ray binaries etc., 
but this group can also include cluster members with in- 



sufficient IR data for an IR classification, either because 
they are highly obscured or due to source confusion in the 
densest part of the cluster. As expected these sources are 
distributed across the field, but a close cluster is seen in 
the center of the cloud (triangles), where Ay > 10 mag. 
Given the low density of background sources in the rest 
of the field, this is statistically a highly significant clus- 
tering. It is very unlikely to find a clustered group of 
unrelated background sources such as AGN exactly be- 
hind the densest part of IRAS 20050+2720. Thus, we 
assume that (most of) these sources must be young stars. 
Only one source in the group is detected in IPHAS, but 
without excess Ha emission. 8 out of 18 sources have no 
IR counterpart at all, 10 are detected in one or more IR 
bands. Mostly the information is insufficient for the clas- 
sification scheme used above. The remaining two sources 
arc classified as stars with contamination by diffuse PAH 
emission in the aperture. They are otherwise compati- 
ble with a class I or II source. All objects in this group 
are brighter for longer wavelengths for all bands where 
they are detected. This is compatible with (reddened) 
YSOs, but their exact status cannot be confirmed with 
the available data. In the following, X-ray sources in this 
region (301.76 < a < 301.81 and 27.47 < 5 < 27.49 - tri- 
angles in Fig. [3]) are marked as X-ray (detected) YSOs = 
XYSO and included in the YSO sample in table [3] and 
[H The cutoff at Ay > 10 mag is chosen to include the 
group of X-ray sources in the center of the cluster. It is 
possible that individual YSOs at lower values of Ay are 
missed, however the only way to identify them without 
sufficient optical or IR information is as a spatial over- 
density compared with the background. 

Hard X-rays, like radio radiation, penetrate much 
higher absorbing column densities than optical or IR pho- 
tons and therefore it is not unusual to detect sources in 
the densest cores of young clusters only in these bands 
(e.g. lHughel 120011: iComito et all [20071: iSchneider et~aT1 
20091 iPravdo et al.ll2009T for Cep A. a massive star form- 
ing region at 730 pc, a distance very similar to IRAS 
20050+2720). 

3.5. Sample completeness 

iFeigelson et al.l (|2005l ) derive a scaling law for the ex- 
pected point source sensitivity of Chandra based on the 
COUP observations. Scaled to the distance of IRAS 
20050+2720 and assuming an average absorbing column 
density of N H = 2 x 10 22 cm" 2 (Sect. [4X2)1 our detec- 
tion limit in the center of the field of view, where the 
PSF is small, should be \ogL x = 29 [erg s -1 ]. We re- 
quire at least 20 photons in a single observation to at- 
tempt spectral fitting, which corresponds to a luminosity 
logix = 30 [erg s -1 ]. 

In the Spitzer observations the mean 90% differential 
completness limit over the observed field is 15.3, 15.2, 
14.6, 13.6 and 8.8 mag for channel 3.6, 4.5, 5.8, 8.0 and 
24/im respectively. To derive the limits, artificial sources 
with different magnitudes are inserted into the observed 
mosaics. The limits give the faintest magnitude that al- 
lowed the pip eline to recover 90% or more of the simu- 
lated sources (|Gutermuth et al.l l2005) . These limits vary 
over the field -of-view and they are on average less sen- 
sitive than in lGutermuth et"aLl (|2009f) (who found 17.1, 
16.8, 14.8, 13.8 and 7.0 mag). The new observations have 
similar exposure times, but the field of view now covers 
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large regions outside the cloud, and thus without back- 
ground extinction, so source confusion is a problem at 3.6 
and 4.5/im. However the density of 24 /im sources out- 
side the cloud is low, thus a lower fraction of simulated 
sources suffers from source confusion and the 90% dif- 
ferential completness limit is at a fainter magnitude. In 
the center of the c luster this work is more sensitive than 
iGutermuth et al.l (|2009T ) because the new data reduction 
significantly enhances the sensitivity in the vicinity of 
bright sources (Sect. |3~5|) . 

The sensitivity of the optical data is very non-uniform 
because of instrumental problems near chip edges and 
because of the wide and non-symmetric PSF, which 
prevents the detection of even moderately bright ob- 
jects close to foreground stars. The faint limit of the 
IPHAS catalog in r' is 21 mag, its bright limit is 
13 mag (|Drew et al.ll20Q5[ ). Thus at the distance of IRAS 
20050+2720 the combination of FLWO photometry and 
IPHAS should be complete to about r' = 20 mag, which 
corresponds to a mass of 0.1, 0.5 and 2.2 for A r 0, 5 
and 8 mag, respectively, for a star of 1 Myr (ISiess etahl 
120001) . 

The detection limits are very different from band to 
band. Also, reddening affects the optical observations 
much more than the IR data. Thus these samples have 
different detection biases and will identify a different frac- 
tion of the true number of YSOs in IRAS 20050+2720. 

3.6. Sample summary 

After removal of foreground objects, the sample of 
YSOs in IRAS 20050+2720 now consists of 276 objects. 
Class I*, I, II and II* objects are classified purely by 
their IR properties. The main limitation for this dataset 
is source confusion in dense regions. All these sources 
may be detected in X-rays. We have identified the X- 
ray bright members of the class III population, which 
we assume to be young, disk-less cluster members, and 
we will use the term "class III" for this group in the re- 
mainder of the article for simplicity. This sample is com- 
plete to about \ogLx — 29 [erg s -1 ]. XYSOs arc defined 
as an overdensity of X-ray sources without IR classifica- 
tion. Thus, by definition, isolated X-ray sources are not 
included. The XYSOs are associated with the densest 
part of the cluster and likely consist of embedded YSOs. 

In the following section we often split the sources into 
three groups called class I, II and III. In this case, class 
I includes class I* and class II includes class II*. 

4. RESULTS 

In this se ction we show ou r resu lts. There is some 
overlap with IGutermuth e"t~a l. (2009), where some of the 
IR data are already presented. However, our work covers 
a larger region and wavelength range. 

4.1. Spatial distribution of sources 

Here we characterize the spatial distribution of the 
sources using a minimum spanning tree analysis. A com- 
plete graph (i.e. each source is directly connected to ev- 
ery other source) is constructed, where the sources of in- 
terest are the "nodes" of the graph. For each connection 
between the sources ( "edges of the graph" ) we calculate 
the angular distance between the sources ( "weight of an 
edge"). A subset of the connections is chosen in a way 
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Fig. 4. — MST for all class II sources, including transition disks. 
The sources (nodes of the graph) are marked with squares. The 
gray scale image shows IRAC 8 fim emission. See the electronic 
edition of the Journal for a color version of this figure. 
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Fig. 5. — Distribution of the branch length for different MSTs. 
See the electronic edition of the Journal for a color version of this 
figure. 

that minimizes the total length but keeps all sources con- 
nected. This is called a minimum spanning tree (MST). 
As an example, Fig. @] shows the MST for all class II 
sources. 

Fig. [5] shows the distribution of the branch length for 
MSTs of class I, II and III objects as well as the MST of 
all YSOs in the cluster. For class I sources most branches 
are relatively short; the distribution for class II sources is 
similar. In contrast, the mean branch length of class III 
objects is much longer. Using an MST has some advan- 
tages to the more commonly employed nearest-neighbor 
(NN) distance and its variants (e.g. the distance to the 
fifth-nearest-neighbor - NN5) because the MST naturally 
results in connected clusters. Although both methods 
perform similarly on round shaped clusters, the MST is 
less susceptible to splitting elongated clusters. This is 
potentially an important drawback of NN methods for 
star formation in the elongated filaments of a molecular 
cloud. 

One important restriction to keep in mind is that all 
these methods analyze the cluster as a two-dimensional 
structure on the plane of the sky. Thus we only see the 
projected separation of sources which is a lower limit to 
the real distance. Also, physically distinct subclusters 
may overlap when projected on the two-dimensional sky. 



In this case mathematical tools such as the MST or the 
NN method cannot provide a distinction between these 
two clusters. 

The branch length distribution for class III objects 
looks qualitatively different from the distribution of 
class I and II sources. The break between two differ- 
ent linear regimes is less pronounced. Also, the average 
branch length is longer than for the other classes, which 
is partly due to the lower number of sources. There are 
70, 188 and 54 class I, II and III sources respectively. 
For a uniform distribution of sources over the same area 
we would expect that the branch length of the class II 
sources is the smallest, and the other classes should have 
distances about a factor of two larger. Instead, class I 
sources have the smallest branch lengths, thus they must 
be most clustered. Class II sources are very similar and 
class III objects have much larger branch lengths, indi- 
cating a more uniform distribution. 

The upper and lower end of the branch length distri- 
bution for class I and II sources can be fitted by two 
lines. The point of intersection for these lines can be 
used to algorithmically define a threshold branch length 
for all stars in a cluster. We cut all edges in the MST, 
which are longer than this threshold branch length of 45", 
leaving several u n conne cted subclusters called "trees". 
IGutermuth et all {2009) introduced this technique in a 
comparison study of young clusters within 1 kpc of the 
sun and they show that it is well suited to a range of 
cluster morphologies. Cutting all branches in the MST 
longer than the threshold branch lengths of 45" (0.15 pc), 
results in a minimum spanning "forest" . The two biggest 
trees have 193 and 37 members, 6 trees have between 3 
and 8 members. Fig.[B]shows the distribution of all YSOs 
in the cluster. Overlaid are the MSTs for the two biggest 
cluster cores ("trees"), which are very different in their 
composition. The bigger cluster core on the western side 
of the cluster ((a,<5) J20 oo = 20 h 07 m 00 s , +27°30', from 
here on cluster core W) has 175 members and a ratio 
of class I:II:III sources of 61:97:17. The smaller, east- 
ern, one ((a,5)j2ooo = 20 h 07 m 40 s , +27°33', from here 
on cluster core E) has 37 members and a ratio of class 
I:II:III sources of 0:31:6, i.e. the ratio of class II to III is 
very similar, but class I sources are absent in the cluster 
core E. This is somewhat surprising given that cluster 
core E is surrounded by a reservoir of gas which can be 
seen in emission in 8/im in Fig. [6] The class I:II:III ra- 
tios reported here are the detected ratios, which do not 
accurately represent the true ratios in the the cluster, 
because class I and II sources are detected from the IR, 
whereas class III objects also require an X-ray detection 
and thus are subject to a different selection bias as ex- 
plained in section 13.51 The number of class III objects 
in cluster core E is potentially underestimated, because 
it is further away from the Chandra aimpoint than clus- 
ter core W. Together with the lower number of class 
sources this argues for a higher age. 

Cluster core E is ou tside of the area surveyed by 
IGutermuth et al.l (|2009f) . but cluster core W is present 
in their analysis. Using the same technique they find a 
lower value for the threshold branch length (28") and 
our cluster core W is divided into two parts, which they 
call Core-1 and Core-2. Even when we cut our MST at 
the same threshold branch length we still only recover a 
single core in the west; Core-1 and Core-2 appear united 
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Fig. 6. — Distribution of all YSOs in the sample. Members of 
the two cluster cores are connected by an MST. The gray scale 
image shows IRAC Sfim emission. See the electronic edition of the 
Journal for a color version of this figure. The black dashed line 
shows one (arbitrary) way to divide the western cluster core in two 
pieces, see text for discussion. 

with the added YSOs reported here. 

However, it is possible to separate cluster core W man- 
ually. This cut may or may not be physical. We present 
it here to discuss limitations of the MST technique. If 
we cut the cluster core W along the dashed black line 
(Fig. [5]), we end up with two subcores. Their ratio of 
class I:II:III is very similar with 21:41:12 and 29:54:7 for 
the northern and southern subcore of cluster core W, re- 
spectively. In the north there are more class III objects, 
but, given the low numbers, this is only marginally sig- 
nificant. Both subcores differ significantly in the mean 
branch length. This is 18" in the northern subcore of 
cluster core W, but only 9.6" in the southern subcore 
of cluster core W. The MST analysis cannot distinguish 
if these are two separate subcores, which just happen 
to overlap in projection on the plane of the sky or if the 
northern subcore of cluster core W is just an extension of 
the southern subcore of cluster core W with lower stellar 
densities. If this is the same physical cluster core, than 
it is asymmetric extending only 2' towards the south of 
the densest core but covering a much larger area in the 
north. 

In summary, we see that class I sources are the most 
clustered and class III objects are least clustered group 
in IRAS 20050+2720. The whole region consists of clus- 
ter core E with 37 members, that is surrounded by gas 
and dust, but has no class I members and the much big- 
ger cluster core W, which contains the densest and most 
extincted parts of IRAS 20050+2720 and has sources of 
all evolutionary stages. 

Table [S] describes the full MST and the MSTs for clus- 
ter cores W and E in order to facilitate a comparison with 
the sampl e of all young stellar clus ters within 1 pc of the 
sun from IGutermuth et al.l (|2009f) . R C irc is the radius 
of the smallest circle that contains all cluster members. 
Because most clusters are not circular we construct a 
convex hull for a cluster and define Rhuii as the radius of 
a circle, that has the same area as the convex hull, with 
an adjustment acc ounting for the number of stars that 
span up the hull (jGutermuth et al"]l2009t ). The aspect 
ratio is R C irc/ Rhuii- In addition, the average density of 
YSOs within the convex polygon a mean is calculated and 
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Fig. 7.— CMD for all YSOs in our sample in SDSS bands. All 
stars are dereddened with the Ajf found during the source clas- 
sification. The plotted errorbars represent only the photometric 
errors in r'i' and do not include any statistic or systematic uncer- 
tainty in Ajf- Ov erplotted are pre-main sequence isochrones from 
ISiess et al.l H2000l '> . See the electronic edition of the Journal for a 
color version of this figure. 

the peak density a pea k o f the cluster after an a daptive 
smoothing (for details see lGutermuth e~a l. 2009). Last, 
the table shows the median nearest-neighboor distance. 

IRA S 20050+2720 was alread y one of the largest clus- 
ters in iGutermuth et al.l (|2009[ ). Here, we have included 
the larger Spitzer field of view and consequently, the size 
of the cluster has increased. Even cluster cores E and 
W separately have a large number of members compared 
to the clusters in the sample. Also, the mean and the 
peak densities ar e signi ficantly higher than reported by 
IGutermuth et al.l (jl009). This is due to the added Chan- 
dra data with allows us to identify cluster members by 
their X-ray emission. Since the distribution of class III 
sources overlaps in space with class I and II sources the 
average density is higher. The XYSOs found in the re- 
gion with the highest extinction are closest together and 
thus define the peak density. This also leads to a very 
small median NN2 distance. 

4.2. Age of the cluster 

While iChen et al.1 (fl997h estimate an average cluster 
age of 1 Myr, they interpret the molecular outflow and 
the stronger reddening in the central parts of subcluster 
W as signatures for a more recent star formation event in 
the last 0.1 Myr. Similarly, the ratio of sources in differ- 
ent stages suggests that cluster core W is younger than 
cluster core E. Further, cluster core W contains many 
class I sources as well as a group of deeply embedded 
objects (XYSOs), both indicative of very young age. 

Fig. [7] shows a r' vs. (r' - i') CMD. It includes all 
YSOs with IPHAS detections and additional nine YSOs 
detected in RI with FLWO. The IPHAS filters are sim- 
ilar, but not identical to, the SDSS filters. We use the 
following transformations which are appropriate for mod- 
erately red stars and have been tested on a G2V star (M. 
Guarcello, private communication): r' = r' IPHAS + 0.12 



L IPHAS 



+ 0.407 + 0.11 xr', 



IPHAS 



. Stars with RI 



and i' 

photometry from FWLO are converted to ri based on the 
inverted relations from Lupton 1 . The transformation be- 
tween ri and r'i' is done according to the formulas given 



on the SDSS website 2 . All stars are dereddened with the 
Ak value, that was estimated during the source identi- 
fication. The ratios of the re ddending coefficients were 
taken from the compilation of Schl egel et al.l (|1998f) . 

The observed scatter in the CMD is, to some extent, 
due to instrumental uncertainties, but also to intrinsic 
scatter in the luminosity caused by bin aries, variability 
and accretion. [Burning ham et al.l (|2005l ) observed a time 
series of photometry for two OB associations and found 
that the temporal variability can account for up to half 
of t he spread in their C MD. In the Upper Sco associa- 
tion [slesnicF^eF^i] (|2008h show with Monte-Carlo simu- 
lations that all their observed scatter can be explained 
by variability, measurement uncertainty, and unresolved 
binaries despite the fact that s pectroscopy is available to 
determine spectral types. Also. iHillenbr and et al.l (2008) 
compared different sets of pre-main sequence evolution- 
ary tracks and found significant offsets between th e ages 
derived from them. Even worse, ICovev et al.l (|2010f ) show 
an example where not even the relative age of different 
star forming regions is conserved when using different 
simulated tracks. They also highlight that most tracks 
do not include the additional luminosity from accretion 
spots, while this contributes to the observed photometry. 
This effect can shift ages by 50-100% or introduce scat- 
ter of the same order, if different objects have varying 
accretion rates. All this will cause an apparent spread in 
luminosity. In addition, even a real spread in intrinsic lu- 
minosity does not require a spread in age. Instead, it can 
be caused by the individual history of st rong accretion 
events for very young stars (Jeffr ies| l2011f) . 

Essentially all stars in Fig. [7] are younger than 3 Myrs 
for cluster core W, possibly with an age spread. Clus- 
ter core E is compatible with an age below 1 Myr, al- 
though no class I source is found there. If there is a 
real age spread in cluster core W, this can explain both 
the existance of ver y young and embedd e d objects like 
jet d r iving sources (iBachiller et all 119951 : ICodella et al.l 
119991: iBeltran et all |2008|), XYSOs^ and class I sources 
and also the older ages indicated in Fig. [71 since the 
younger members are likly too embedded to be observed 
in r'i' and thus would not be shown in the figure. An 
alternative explanation could be the presence of a mech- 
anism that destroys disks, such as radiation and winds 
from a massive star. However, no massive star is ob- 
served in this region; also, cluster core E is surrounded 
by much stronger gas emission than cluster core W, so 
any gas destruction would have to be confined to the 
disks, but not the surrounding material. 

No systematic differences are seen between class I and 
II objects, if they are plotted separately in the figure. 

We can confirm a cluster age of < 3 Myrs for 
IRAS 20050+2720 and < 1 Myr for cluster core E, al- 
though the errorbars in Fig. [7] do not contain all sources 
of uncertainty. While we cannot distinguish between a 
past star formation event and ongoing star formation 
with certainty, the scatter in luminosity and thus ap- 
parent age is large and could indicate a real age spread 
in cluster core W. 

4.3. X-ray properties 



1 http:/ /www. sdss.org/dr7/algorithms/ 
sdssUBVRITransform.html#Lupton2005 



2 http:/ /www. sdss.org/dr7/algorithms/ 
jeg_photometric_eq_drl.html#usno2SDSS 
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TABLE 5 

Geometrical properties of IRAS 20050+2720 and its cluster cores. 
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Fig. 8. — Histogram of dereddened K band magnitude for class II 
sources. 



4.3.1. Detection fraction 

Table 0] shows the number of IR sources and X-ray 
sources in each class. Few class I sources have X-ray 
detections, but for class II sources we can compare the 
properties of those with and without an X-ray detec- 
tion. Figure |S] shows their K band magnitude distribu- 
tion. All values are dereddened with Ak as found during 
the source classification. The mean value is Ak = 0.6 
over all class II sources. X-ray bright sources are also 
bright in the 2MASS and Spitzer/IRAG bands. Half of 
all sources with tuk < 12 mag are X-ray sources but 
the detection fraction drops sharply for less luminous 
sources. In a magnitude limited sample with dered- 
dened rtiK < 13 mag the Ak values are fully com- 
patible between X-ray detected and X-ray non-detected 
sources (KS-test probability that Ak is drawn from the 
same sample p — 0.99), but the magnitudes are not 
(p = 1~ 10 ). This indicates that the non-detection of 
X-rays is not caused by higher extinction but is intrinsic 
to the sources, i.e. objects fainter in the IR are weaker X- 
ray emitters and thus remain undetected. A ccording to 
pre-m ain sequence evolutionary tracks from ISiess et al.l 
(2000) niK = 12 mag corresponds to stellar masses of 
0.3, 0.7 and 0.9 M Q at 1, 2 and 3 Myrs, respectively. 

The regions of the highest absorbing column densities 
are all part of cluster core W, thus this should show the 
lowest X-ray detection fraction, if the luminosities are 
uniform, contrary to the observed distribution. Also, 
stars in cluster core W are brighter in K band than 
sources in the other regions. This might indicate an 
evolutionary trend: Cluster core W still contains many 
class I sources and also the class II sources could on aver- 
age be younger and less contracted. This would explain 
the higher K band luminosity a nd also the higher X-ray 
flux, since Wi nston et al.l (|2010t ) show that the X-ray flux 




2 3 4 

N u [10 22 cm- 2 ] 

Fig. 9. — Histogram of absorbing column densities Njj for class I, 
II and III sources. See the electronic edition of the Journal for a 
color version of this figure. 

is linearly related to the surface area, such that larger 
YSOs have a higher luminosity. 

4.3.2. Absorbing column density 

Some cluster members have been identified by their 
X-ray emission in Sect. I3.4[ but for an IR selected sam- 
ple table @] shows that the X-ray detection fraction in- 
creases with the evolutionary status of the sources from 
class I* to class II*. This is expected because the enve- 
lope of younger objects absorbs the soft X-rays and thus 
reduces the observable flux. In the following, we present 
results for fits to the X-ray spectra (see Sect. I2.2|) . Be- 
cause of the low signal the fit does not converge for all 
sources and we ignore those fits, where the statistical 
error on the fit value includes the complete parameter 
space. Fig. [9] confirms that class I sources are indeed 
more embedded than class II and III sources. According 
to a KS test, the probability that class I and II sources 
are drawn from the same sample is 23%. The probability 
that class II and III sources are part of the same sample 
is only 3%. This has also be en observed in other s tar 
foming regi ons e.g. the ONC dP rismzano "eTal1[2008l) or 
NGC 1333 (jWinston et al.ll2010t) . 

While the deeply embedded sources may remain unde- 
tected, so that the figure does not represent the complete 
distribution of extinctions, the differences in the distri- 
bution of Nh between different classes are real. Both 
class I objects with fits have Nh — 2 — 6 x 10 22 cm -2 , 
much greater than the average class II or III source. Less 
absorbed sources are easier to detect and if less absorbed 
class I sources were present, they would have been de- 
tected. The lower detection fraction of younger sources 
in X-rays (table 0J means, that class I sources are fainter 
or they are so deeply embedded that the detections in 
Fig. IH1 represent only the less absorbed end of the dis- 
tribution. Thus, we conclude that sources with larger 



11 




Class I 
Class II 
Class III 
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Fig. 10. — Histogram of X-ray temperatures kT for class I, II 
and III sources. Sources with jVjy > 10 22 cm -2 are shown in 
hatched histograms. No separate histogram is plotted for class I 
sources, because all of them show Nh > 10 22 cm~ 2 anyway. See 
the electronic edition of the Journal for a color version of this figure. 

infrared excesses are indeed more absorbed in X-rays. 

4.3.3. Temperature 

Figure [10] shows the distribution of temperatures for 
the different classes (open histograms). Clearly, class I 
sources are on average hotter than the other sources. 
However, this could be due to a detection bias: Class I 
sources are generally more absorbed, thus soft X-ray 
emission remains undetected. To avoid this problem a 
subsample with similar Nh can be used to compare the 
temperatures, because all sources suffer the same bias 
towards higher temperatures. Unfortunately, the low 
number of sources requires a relatively broad range of 
2 x 10 22 cm -2 < Nh in order to retain a sufficient number 
of sources (hatched histograms). This interval contains 
both class I sources. Although a KS test is not strictly 
valid for such small numbers, it indicates that the dis- 
tribution of the plasma temperatures for class II and III 
sources in this subsample differs from the distribution 
in the full sample, because lower temperature plasma 
remains undetected. The subsample contains 2 class I 
sources, 6 class II sources and 3 class III sources. Even 
in this small subsample a KS test shows that the three 
temperature distributions are incompatible. The class I 
sources are hotter than the class II sources and those in 
turn are hotter than the class III soures. However, there 
may still be a bias, because class II and III sources are 
mostly on the lower end of the Nh range in the subsam- 
ple. Due to the larger source numbers we can repeat this 
test for class II and III sources in smaller Nh bins. For 
N H = - 1 x 10 22 cm- 2 , N H = 1 - 2 x 10 22 cm" 2 and 
Nh = 2 — 3 x 10 22 cm -2 the mean temperature of the 
class II sources is always higher than for class III sources 
and the probabilities that the temperatures for the two 
classes are drawn for the same distribution are <1%, 3% 
and 15%, respectively. 

Thus, class II sources in our sample seem to be hotter 
on average than class III sources and it seems likely that 
class I sources are even hotter than class II sources. This 
has been seen before in other young cl usters. While no 
signifi cant trend is found in the ONC iPrisinzano et al.l 
(2008) show in their Fig. 13 that class II objects are 
margi nally hotter than class III sources. IWinston et al.l 
(2010) studied Serpens and NGC 1333. In these two 



clusters, the class III sources are found to have lower 
values of kT than the class II. Both samples are thought 
to have similar completeness for each class. The K-S 
probabilities that the class II and III temperature arise 
from the same distribution are 21% (Serpens) and 3.3% 
(NGC 1333). Our result agree with these findings, and 
the fact that we can show the effect in different Nh 
ranges makes it very unlikely that this is caused by ab- 
sorption only. However, the spectral models we use are 
based on a single temperature. This is an adequate de- 
scription of the observed data with its low signal, but we 
know from observations of closer CTTS and WTTS that 
the real temperature structure in coronae is more com- 
plex. Also, the comparison in different Nh bins could 
be flawed, if the Nh is caused by a partial absorber in 
class II or III sources, because only a part of the X-ray 
emission is covered by accretion funnels. Both consider- 
ations cannot be tested with the given data quality, thus 
we caution that the observed difference in temperature 
might not represent the real coronal temperature very 
well. 

On first sight, the observed temperature trend seems 
to contradict the soft exces s found in CTTS, mostly 
from O VII line emission dRobrade fc Schmittl [2007; 
iGudel fc Telieschl [20071: iGuntherl 12011ft . but Chan- 
dra/AClS is less sensitive at low energies than XMM- 
Newton and only a very small absorbing column density 
(Nh « lx 10 21 cm -2 ) is needed to absorb the soft emis- 
sion. Thus, the data in Fig. [TOlshow that the hot plasma 
component is hotter in class I and II sources than in 
class III sources, but the data do not constrain the pres- 
ence of a soft emission component from accretion shocks. 

4.3.4. Luminosity 

In clusters of stars there is a wide distribution of X-ray 
luminosities. Luminosities for stars with spectral fits are 
given in table [2j but this applies only to a fraction of the 
YSOs with X-ray detections. To estimate a rough lumi- 
nosity for the remaining sources, we fit a relation between 
count rate c and computed unabsorbed source flux / to 
the YSOs with spectral fits. We find the following rela- 
tion: / = c* (4±2) x 10 -9 erg s _1 cm -2 count -1 . While 
this is very uncertain for a single star because the pro- 
portionality constant really depends on the temperature 
and absorption, this fit is a reasonable approximation for 
a sample of stars, because the relation represents the av- 
erage values of Nh and kT for those YSOs with spectral 
fits. Figure [5] indicates that the X-ray detected sources, 
at least for class II objects, are the brighter in K-band 
and thus are likely more massive stars. In this case, the 
sources with weak X-ray emission might just be the less 
massive and thus intrinsically fainter than their counter- 
parts with X-ray spectral models, but could still have 
similar Nh and kT values. 

Above \ogLx = 29.8 the X-ray luminosity function 
(XLF) is dominated by sources with spectral fits, thus 
the sample is fairly complete and we are confident in 
the Lx values. However, soft sources could be hidden 
in the dense parts of the cloud, even if they have a 
high Lx- Most sources below this value have luminosi- 
ties which are only estimated from the count rate. Fig- 
ure [TT] shows the cumulative X-ray luminosity function 
for IRAS 20050+2720 and several other star forming re- 
gions for comparison: The Orion nebula cluster (ONC, 
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Fig. 11. — The plot shows cumulative X-ray luminosity func- 
tions (XLF) for different star forming regions normalized to 1.0 
at logix = 29.3. The data for Serpens and NGC 1333, which 
have less than 100 measured values, are shown as individual points. 
The shaded area shows an analytical lognormal distri bution with 
fi = 29.3 and <r = 1.0 as suggested by Fcigclson ct al. (2005). See 
text for data sources. See the electronic edition of the Journal for 
a color version of this figure. 



Getman et all l2005f ). IC 348 (jPreibisch fc Zinneckerl 
2001 . and Serpen s and N GC 1333 (jWinston et al]l2lH0h . 

iFeigelson et al] (|2005l ) suggested that the XLF of 
young star forming regions could be universally described 
by a lognormal distribution with /i = 29.3 and rr = 1.0. 
In Fig. 1111 The remaining clusters in the figure are com- 
plete only for higher luminosities. The high luminosity 
end of the distribution is defined by very few massive 
O and B stars, which are not present in the lower mass 
clusters. The ONC follows a lognormal curve very well 
here. Differences on the high-mass end of the XLF are 
seen betwee n the ONC and NG C 2244, a similarly mas- 
sive cluster (|Wang et al.ll2008l ). and there seems to be 
an excess of sources around logLx = 29.3 in Cep B 
(jGetman et al.1 120061 ) . Here we compare it to several 
clusters of lower mass. The ONC sample should be 
nearly complete down to M* < 0.2M Q . Serpens and 
NGC 1333 have the steepest XLF and the lowest cluster 
mass (< 200 M m ), IC 348 and IRAS 20050+2720 (300 
and 430 M Q , iGutermuth et all I2009T) have a shallower 
XLF and intermediate masses and the ONC, the most 
massive cluster in the figure, has the shallowest slope 
and follows a lognormal distribution with a = 1.0 most 
closely. We compare the XLF above \ogLx — 29.8 with 
2-sided KS tests. In Serpens there are only 10 sources 
in this range, thus it is compatible with any other distri- 
bution. Apart from Serpens, the chance that the ONC 
sample and any other XLF are drawn from the same 
parent distribution is < 2%. IRAS 20050+2720 agrees 
with IC 348 and NGC 1333 on the 30% level, while the 
chance that IC 348 and NGC 1333 are drawn from the 
same parent sample is 19%. Thus, we see a devitation of 
the observed XLFs for low-mass clusters from the COUP 
XLF, a further indication that the XLF is not universal. 

The more massive clusters contain a larger number of 
massive stars of the spectral types B to F, which are 
by themselves relatively faint in X-rays. However, a 
high fraction of them are binaries with an X-ray bright 
companion, i.e. a G or K type star. These binaries 
could provide additional X-ray luminosity in the range 
log Lx = 29 — 30 and thus tilt the observed distribution 
to explain why more massive clusters have a shallower 
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Fig. 12. — Relation between gas column density Nh and ex- 
tinction Ak- Different symbols show differ ent source types for 
our sa mple. For comparison, all CTTS from [Gunthcr & Schmittl 
are also shown. Lines indicate a fit to the data for IRAS 
20050 +2720, the interstellar gas-to-dust r atio from IVuong et al] 
l(2003T l and the ratio found in NGC 1333 bv lWinston et al.l l|20fc 



I). 

All lines have the same slope in this log-log plot, because they are 
all derived assuming Njj oc Aj^ and only differ in the proportion- 
ality constant. See the electronic edition of the Journal for a color 
version of this figure. 

XLF. 

4.4. Gas absorption and dust extinction 

The absorption of X-rays is mainly due to metal atoms 
and ions in the line of sight. It is generally expressed as 
the equivalent hydrogen column density Nh- In con- 
trast, optical and infrared extinction is mainly caused by 
dust. Simplistically, it should be possible to calculate the 
gas-to-dust ratio along the line of sight for those sources 
with an X-ray fit and sufficient IR information to esti- 
mate Ak- However, the situation is complicated by grain 
coagulation and ice growth. 

Figure [H shows the Nh vs. A k for IRAS 20050+2720 
with different symbols for different evolutionary classes. 
Due to the uncertainty in the intrinsic color of the stars, 
there is an uncertainty of 0.1 in Ak (R- Gutermuth, 
personal communication). A linear fit gives the rela- 
tion Nh = (1-1 ± 0.1) x 10 22 A K cm" 2 mag" 1 . There is 
no significant difference between class II and III sources 
and separate fits for cluster core E and W lead to com- 
patible values. The Nh/Ak ratio is lower than the 
interstellar ratio of Nh = 1-8 x 10 22 Ak cm" 2 mag -1 
(|Predehl fc Schmittl Il995t ). where A v was transformed 
to Ak with the tabulated Ak/Av of ISchlegel et al.1 
(1998). However, values lower than in IRAS 20050+2720 



are sometimes seen in young star forming regions, e.g. 
N H = (0.89 ± 0.13) x l O 22 ^ cm" 2 mag" 1 for class III 
sources in NGC 1333 ([Winston et al.ll2010f ) and N H = 
(0.63 + 0.23) x 10 22 A K cm" 2 m ag" 1 for class III so urces 
in Serpens. In the p Oph cloud IVuong et all (|2003h find 
1.4 + 0.1 x 10 22 ^a- cm" 2 mag" 1 . 

The situation appears more consistent with the ISM 
in higher mass star formi ng regions. In RCW 108 
(with at least one 08 star) IWolk et al.1 (2008) obtained 
N H = (1-6 - 2.60) x lO 22 ^ cm" 2 mag" 1 for differ- 
ent subregions of the cl uster. For the 05 bearing star 
forming region RCW 38 IWolk et all (l200l find N H = 
1.8 x 10 22 Ak cm -2 mag" 1 , again in agreement with the 
interstellar ratio. 

In an effort to understand this bifurcation, figure [12] 
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also shows Nu vs. Ak for the very different regime of 
close- by and little absorbed class I I sources from the sam- 
ple of iGiinther fc Schmittl (2008, see references therein 
for data sources). Most of those sources have been ob- 
served with X-ray grating spectroscopy, thus the Nh val- 
ues are more reliable even for small column densities. On 
the other hand, the Ay values for this sample have been 
collected from the literature and are derived with differ- 
ent methods. In the plot sources with a literature value 
of Ay — are shown at Ak = 0.01 mag, because it 
is difficult to measure a small reddening. Most of these 
sources are found significantly above the Nh/Ak of the 
d iffuse ISM. 

IWinston et al.l (|2007l ) proposed three possible reasons 
for the low Nh/Ak values. First, the X-ray temper- 
atures could be systematically underestimated resulting 
in an over estimate of the Nh values. Only X-ray grating 
spectroscopy can break the Nh to kT degeneracy in the 
fits. However, this is not feasible for class II or III sources 
beyond about 150 pc. Second, that the value of Ak has 
been overestimated due to changes in the form of the red- 
dening law; however, the required change would be much 
greater than has so far been observed. The third reason 
was that the growth and coagulation of grains would in- 
crease the extinction per unit grain mass in the cloud. 
They find the first two reasons unlikely, and argue that 
grain growth is the most likely explanation. 

While both Nh and Ak are measures of the mate- 
rial along the line of sight, Nh is essentially a mea- 
sure of number of atoms while extinction is affected by 
particle shape and the grain size distribution. In the 
case of a star, which is surrounded by a disk and has 
a gas accretion column, there are various layers which 
can have different particle shapes and grain size distri- 
butions. Nonetheless, in IRAS 20050+2720 class II and 
III sources show the same Nh I Ak ratio, although class II 
sources have more circumstellar material. Thus, either 
the extinction laws of the circumstellar material and the 
cloud are similar or the circumstellar extinction is not a 
significant contributor to the total extinction. 

If we view Nh/Ak as a gas-to-dust ratio, IRAS 
20050+2720, Serpens and NGC 1333 are dust rich (or gas 
poor), compared with the interstellar extin c tion. In con- 
trast, the sample from IGiinther fe Schmittl (|2008| ) shows 
most sources to be dust depleted (or gas enriched) with 
much lower absolute values for the extinction. These 
sources are mostly from the Taurus-Auriga molecular 
cloud complex, which does not have the high density 
cluster cores seen in IRAS 20050+2720, Serpens and 
NGC 1333, and contributes little to the line-of-sight ex- 
tinction. Thus, their extinction should be dominated by 
circumstellar material, which could be dust-depleted due 
to the energy input from the star, gravitational settling 
of grains in the disk-midplane and planet formation. 

In contrast, in the younger star forming regions, we 
mostly see the cloud material. Although we cannot rule 
out that these clouds could be dust-rich (perhaps from 
a high abundance of carbon) or poor in X-ray absorbing 
metals, there is no independent evidence for these expla- 
nations. Alternatively, grain growth and coagulation can 
lead to an increase in the extinct i on per Hydrogen atom 
in th e near IR-bands (lJural fl980: We ingartner fe Draind 
1200 ll) . IChiar etlrtl (|2007l) report a decrease in the ratio 
of the depth of the 9.7 /im silicate feature to the color 



excess E(J — K) in molecular clouds relative to that in 
the diffuse ISM. They note that models predict grain 
growth can increase the amount of extinction per Hy- 
drogen atom in the near-IR wit hout strongly affecting 
the s hape of the reddening curve (|Ossenkopf fc Henningl 
1994); note that the value of Ak is calculated primarily 
with the E(J — H) color excess and E(H — K) color ex- 
cesses; these should follow a similar trend as the E(J—K) 
color excess. Models of the grain coagulation in molecu- 
lar clouds bv lOrmel et al.l ()201lD show the E(J — K) per 
unit column density can increase relative to the diffuse 
ISM due to coagulation until the grains start exceeding 
1 pm in size. 

In this scenario, the low Nh/Ak ratio is due to 
grain growth in the cold, dense molecular cloud mate- 
rial surrounding low-mass star forming regions; although 
there is some dispersion in the amount of grain growth 
since IRAS 20050+2720 shows a higher N H /A K than 
NGC 1333 or Serpens. In contrast, the extinctions and 
gas column densities for class II objects in the Taurus- 
Auriga region may have a significant component from the 
circumstellar disk; grains in these disks may have grown 
to sizes in the mm range, too large to cause s ubstantial 
reddening, thus the IGiinther fc Schmittl (|2008[ ) measure- 
ments show a larger Nh/Ak ratio. Alternatively, scat- 
tered light from the grains may also contribute to the 
colors of the objects causing an underestimate of the ex- 
tinction. The RCW 108 and RCW 38 are massive star 
forming regions, whose members could supply sufficient 
radiant flux and/or mechanical power to evaporate grain 
mantles an d break up coagu lated dust in large parts of 
the cloud (|Smith et al.lll999D . leading to an extinction 
law closer to th e interstellar one. 

Alternatively, IVuong et al.l (|2003f) show that the low 
Nh/Ak ratio in p Oph, which is compatible with the 
value found in IRAS20050+2720, can be described by a 
change in the local abundance, that is consistent with 
recent revisions of the local (solar) abundance (see their 
paper for a detailed comparison of different sets of abun- 
dances). However, this would require a different set of 
abundances for each of p Oph and IRAS 20050+2720, 
Serpens, and NGC 1333, while RCW 108 and RCW 38 
both should have the same abundances as the galactic 
ISM. The abundances in the stellar photosphere have 
not been measured to the precision required to compare 
the different clusters, thus this idea cannot be confirmed 
with the available data. 

In summary, the cold molecular clouds of low-mass star 
forming regions have Nh/Ak ratios that are systemati- 
cally lower than that of the diffuse ISM due to their cold, 
dense (> 10 4 cm -3 ) and relatively unperturbed environ- 
ments which facilitates the growth of grains. 

5. SUMMARY 

We present multi-wavelength data for the young star 
forming region IRAS 20050+2720 from the Spitzer and 
Chandra satellites and the ground-based 2MASS and 
IPHAS surveys as well as targeted FLWO photometry. 
YSOs of class I and II are identified based on their IR 
properties. X-ray sources with optical or IR counterparts 
that are compatible with the cluster distance are added 
to this sample as well as a group of X-ray sources that 
are identified as cluster members based on their close 
association with the densest cluster core. 
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We use a minimum spanning tree (MST) analysis to 
characterize the spatial distribution of sources in IRAS 
20050+2720. Overall, class I sources are more clustered 
than class II sources, and class II sources are more clus- 
tered than class III sources. 

Two main cluster cores can be identified in IRAS 
20050+2720. Cluster core E consists only of class II 
and III sources, but seems younger in the r' vs. r' — i! 
diagram. It is surrounded by an IR nebulosity. Clus- 
ter core W has a dense and highly extincted core in the 
South and a much less dense region in the north. It is 
not clear if those two regions are distinct groupings along 
overlapping lines-of-sight, or if they are part of the same 
physical cluster core. The age of IRAS 20050+2720 is 
< 3 Myrs and it might contain populations of different 
ages. 

In X-rays we detect preferentially the more massive 
YSOs; the detection fraction is higher for more evolved 
evolutionary stages, where the absorbing column densi- 
ties are lower. Still, the detected class I objects are both 
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